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Non-Gaussian statistics of electrostatic fluctuations of hydration shells 
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We report the statistics of electric field fluctuations produced by SPC/E water inside a Kihara solute given 
as a hard-sphere core with a Lennard- Jones layer at its surface. The statistics of electric field fluctuations, 
obtained from numerical simulations, are studied as a function of the magnitude of a point dipole placed 
close to the solute-water interface. The free energy surface as a function of the electric field projected on the 
dipole direction shows a cross-over with the increasing dipole magnitude. While it is a single-well harmonic 
function at low dipole values, it becomes a double-well surface at intermediate dipole moment magnitudes, 
transforming to a single-well surface, with a non-zero minimum position, at still higher dipoles. A broad 
intermediate region where the interfacial waters fluctuate between the two minima is characterized by intense 
field fluctuations, with non-Gaussian statistics and the variance far exceeding the linear-response expectations. 
The excited state of the surface water is found to be lifted above the ground state by the energy required 
to break approximately two hydrogen bonds. This state is pulled down in energy by the external electric 
field of the solute dipole, making it readily accessible to thermal excitations. The excited state is a localized 
surface defect in the hydrogen-bond network creating a stress in the nearby network, but otherwise relatively 
localized in the region closest to the solute dipole. 

Keywords: Hydration, surface phase transition, hydrogen bonds, optical spectroscopy, solvent effect 



I. INTRODUCTION 

Polarization of polar liquids by external electric fields 
is one of classical problems of condensed matter theoryi^ 
The ideas advanced by Debye and Langevine assign the 
average dipole moment in a liquid exposed to an ex- 
ternal electric field to a function increasing monoton- 
ically and linearly (Debye) or non-linearly, with satu- 
ration (Langevine) with the increasing field strengthi^ 
The polarized medium, in return, produces a field of 
its own which, combined with the external field, yields 
the macroscopic Maxwell field. One thus arrives at 
the concept of dielectric screening implying the ability, 
quantified by the macroscopic dielectric constant, of the 
medium to polarize and produce its own electric field, in 
response to an external perturbation. 

The electric field of a polar liquid can be alternatively 
probed by an external multipole. The electric field R(r) 
inside a liquid is a fiuctuating variable depending on the 
instantaneous liquid configuration represented by a point 
in its phase space F. Thus an external dipole moment mo 
placed outside or inside the liquid will gain the instanta- 
neous electric energy e — — mg • R. The liquid, in turn, 
is polarized by mg such that its field at the position of 
the dipole fluctuates around the average field Rq, which 
is aligned along mo and is known as the Onsager reac- 
tion field,i This average field is accessible experimentally 
from solvent-induced shifts of optical dyesi^i^ 

A fluctuation of the electric fleld SR{T) = R(r) - 
Ro out of equilibrium requires reversible, non-expansion 
work applied to the systemji this work is quantified by 
the Landau functional -F(R). We will be interested here 



only in the projection of the field on the direction iho = 
mo/mo of the external dipole and thus will set up the 
"reaction coordinate" or "order parameter" R = iho ■ R.— 
The free energy surface can then be found by standard 
prescriptions^ singling out the order parameter from the 
manifold of the system degrees of freedom 



5{R~ R{T))e-'^"dr. 



(1) 



In this equation, H is the Hamiltonian of the solution 
involving a solute carrying no dipole (mo = 0), and 
/3 = l/(k^T) is the inverse temperature. The free energy 
F{R) thus describes thermal fluctuations of the electric 
field in the absence of the dipole mo. 

Empirical observations suggest that the Onsager reac- 
tion field i?o is a linear function of the dipole moment 
inducing it.— This observation can be mathematically 
cast into the requirement of a harmonic form of F{R) 
(Fig. nil) 



F{R) = R^/{2k). 



(2) 



The free energy changes when the dipole moment is 
turned on into 



T{mQ,R) = -moR + F{R). 



(3) 
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The minimization of this free energy in respect to the 
fleld R yields the reaction fleld linear in the solute 
dipole, Rq = KrriQ. This is the familiar linear response 
approximation)^ It articulates two physically signiflcant 
concepts: (i) the response of the medium (here electric 
fleld) grows linearly with an external perturbation (solute 
dipole moment) and (ii) the perturbation does not alter 
the spectrum of the medium fluctuations which therefore 
can be calculated or measured from the properties of the 
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system in the absence of the perturbation (zcro-dipole 
solute in our case). 

The distribution of the solvent electric field P{R) oc 
eKjp[— (3 F (R)] is a Gaussian function with the variance 
cr| = (((5i?)2) equal to k/P, SR = R - Rq. Since fluc- 
tuations of the electric field introduce fluctuations in the 
optical transition energies,—"— the stiffness parameter k 
(also known as response function) is experimentally ac- 
cessible from the inhomogeneous width of an optical line. 
At the same time, since the average energy, Cq = — ktoq, 
and the distribution width are given in terms of the same 
parameter k, one gets the fluctuation-dissipation^ rela- 
tion between the average energy and the energy variance, 
eo = — /?(((5e)^). The parameter 

XG = -l3{{6ef)/eo (4) 

can therefore be considered as the non-linearity parame- 
ter quantifying deviations from the Gaussian statistics of 
the electric field fluctuations xg = 1 for the Gaussian 
statistics when linear response holds and the stiffness pa- 
rameter can be calculated either from the reaction field 
or from the field variance 

K ^ Ro/mo ^ Paji- (5) 

The inverse stiffness parameter has the meaning of 
the characteristic volume in which a small probe dipole 
will be influenced by the fluctuations of the electric field 
R. One can then introduce the characteristic length to 
which the field fluctuations penetrate inside the solute 

iR = i^)-'/' (6) 

This length will then define the thickness of the surface 
layer within which a dipolar probe (e.g., optical dye) will 
be affected by the solvent (i.e., the interaction energy will 
be greater than ksT). 

The harmonic free energy profile [Eq. ([2])] sketched in 
Fig. [IJi is not the only concevable form of the free energy 
surface. The reaction field Rq building up in response to 
a growing dipole moment niQ is produced by the solvent 
dipoles increasingly aligned with the external field. This 
externally forced alignment leads to frustration of the 
local orientational structure within the liquid, which can 
release itself into a new structure represented by a second 
minimum of the Landau functional F{R) (Fig. [Tla). 

Several fundamental consequences might be antici- 
pated if the situation sketched in Fig. [T]d is realized. 
First, the curvature of the free energy at the higher- 
energy minimum is not necessarily equal to the curvature 
at the lower-energy minimum. In other words, the fluctu- 
ation spectrum around the minimum i?i might be altered 
compared to the spectrum in the absence of the external 
field. If this happens, the assumption (ii) listed above (no 
alteration to the fluctuation spectrum) does not apply 
anymore. Further, the free energy Fq at the second min- 
imum should be lifted compared to F(0) = in expec- 
tation of no spontaneous field in the ground state. How- 
ever, the field of the dipole mo will bring the minimum 
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FIG. 1. Sketch of the standard expectation (a) and the pic- 
ture proposed in this study (b) for the free energy of creat- 
ing field _R at a point within the solute when no dipole is 
present [Eq. In our simulations, a point dipole is placed 
at that point, which modifies the free energy to T{mo,R) 
[Eq. ([3])]. The traditional linear-response approximation an- 
ticipates F{R) to be harmonic in the whole range of i?-values 
of interest. In contrast, we suggest the existence of an exci- 
tation state, with the corresponding minimum lifted by the 
free energy Fq. This excited state, too high to be observed in 
the homogeneous liquid or in non-polar solution, can become 
observable by placing the dipole mo into the solute, which 
lowers the free energy gap between the lower and higher min- 
ima by —moRi- The curvature of the excited state near its 
minimum might be different from the corresponding curvature 
at the lower minimum, thus altering the basic assumption of 
the linear-response approximation, the invariance of the spec- 
trum of the solvent fluctuations with respect to the magnitude 
of mo. The overall shift of F{R) in the external field is not 
shown in the diagram; the positions of the lower minimum at 
zero and nonzero dipoles mo are therefore made coincide. 



of F{mo,R) down by — moi?i (Fig. [T]). This external- 
field effect offers a unique opportunity to study excitation 
states of a liquid interface not significantly populated by 
thermal excitations due to a high value of Fq. 

Interfaces of water with solutes of nanometer dimen- 
sion is a particularly suitable candidate for observing sur- 
face excited states. Large solutes break the network of 
hydrogen bonds of bulk water producing unsatisfied sur- 
face bonds. The result is a specific orientational struc- 
ture of the surface waters with their dipoles preferen- 
tially oriented parallel to the interface J^"— While strong 
hydrogen-bond interactions provide a sufficient driving 
force for in-plane orientations, the electric field of exter- 
nal multipoles will compete with this alignment, possibly 
switching it, in a discontinuous fashion, into a new inter- 
facial structure. In fact. Sharp and co-workers have found 
that water at the solute interface occupies two states of 
local structure characterized by straight, "ice-like" hy- 
drogen bonds around non-polar solutes and bent bonds, 
required by the dipoles alignment along the electrostatic 
field of polar/ionic solutes J^ii^ The distinction between 
polar and non-polar solvation can be attributed in this 
picture to changes in occupations of these two hydrogen 
bonding states. What has been still missing from that 
picture is the ability of surface waters to switch between 
the two states via a concerted cross-over. The existence 
of such a cross-over, occurring with increasing electro- 
static field of the solute, is the main finding of this study. 

The system we study here is a non-polar solute im- 
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FIG. 2. Cartoon of the solute-solvent configuration. The so- 
lute is a hard core of radius rns covered with the LJ layer of 
the width aos [Eq. ([TJ]. The distance ros = '"hs +o"0s approxi- 
mately corresponds to the first peak of the solute-solvent pair 
distribution function gos{r) [also see Fig. [S] below]. The point 
dipole mo is placed at the distance Vd from the solute center 
and, correspondingly, the distance r^s = ros — rns from the 
solute-solvent interface. The solvent is SPC/E water at 1 atm 
and 273 K (the melting temperature of SPC/E water is 215 
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FIG. 3. Reaction field i?o (in V/A) vs solute dipole mo for 
RD9 {jd = 9 A) and RDIO (r^ = 10 A) configurations (see 
Fig. [2ll . The dashed lines show the slopes of the linear por- 
tions of Ro{mo). The scaling of Rq with mo is approximately 
quadratic in the intermediate region. The vertical arrow in 
the lower panel indicates the equilibrium field Ri in the ex- 
cited state obtained by extrapolating the reaction field in the 
excited configuration to zero solute dipole (Fig. [1]). 



mersed in SPC/E water— and interacting with it by a 
Kihara potential characterized by a hard-sphere core, 
ths = 9 a, and two Lennard- Jones (LJ) parameters, 
(Tos = 3 A and eo, = 0.65 kJ/mol, 



r - rns 



Cos 



(7) 



We then modify the system by introducing a dipole mo- 
ment at the distance rd from the solute center (Fig. ^ 
and watch the changes in the distribution of the water 
field R experienced by the dipole when the dipole mag- 
nitude grows. While the system follows the standard 
expectations of linear response at low solute dipoles mg, 
the situation changes at some critical value of the solute 
dipole when a new, non-linear, functional form for the re- 
action field emerges (Fig. [3]). The reaction field returns 
to a linear trend with mg when the dipole keeps growing, 
but extrapolates to a non-zero value at mg — >■ 0, as would 
be expected if the system has landed in a high free-energy 
minimum, as in Fig. [T]d. Further, the curvature around 
this second minimum (slope of Rq vs mg in Fig. [3]) also 
experiences alteration, as we have described above. The 
basic physical picture is the emergence of a finite domain 
of surface waters in an excited state different both from 
the bulk and from the rest of the interface. The details of 
the simulation protocol can be found in the Supplemen- 
tary Material (SM}2i and we proceed to the discussion 
of the results. 



II. THREE-STATE MODEL 

Before proceeding to a more detailed description of the 
simulation results, we briefly outline the analytical frame- 



work used here to understand the outcome of the simula- 
tions. This is a phenomenological three-state model rep- 
resenting the simulation data in terms of three distinct 
states of surface waters. 

The model includes three states described by Landau 
free energy functions Fa{R) [Eq. ([T])] with a = g,l,2 
standing for the ground state (g) and two exited states 
corresponding to the minima R^ at the positive (1) and 
negative (2) reaction field values. Two excited states, 
instead of one as in Fig. [l] are introduced to satisfy the 
condition of vanishing reaction field at toq ~^ 0. The 
minima of these states are lifted by the energies Fa from 
the free energy minimum Fg (0) of the ground state. Each 
of Fa (R) curves is modified by the presence of an external 
dipole mg into J-Q(mg, R) — Fa{R) — mgi?. as in Eq. ([3]). 
The global free energy of the system is then found as the 
statistical trace over the three states as follows 



-/3J^(mo,-R) 



-/3J^„(mo,_R) 



(8) 



From MD simulations, we have found little asymme- 
try in the distribution of water's dipolar directions at 
the interface, as judged from the first-order orientational 
parameter 



Pi 



j=l,N' 



(9) 



Here, the sum runs over iV^ waters in the solute's first 
hydration layer, ih-,- = nij/m, and Vj = Yj/rj is the unit 
vector defining the position of the first-shell water rel- 
ative to the solute's center. The order parameter p( is 
close to zero indicating that dipoles oriented parallel and 
antiparallel to the surface normal occur with approxi- 
mately equal probability. Therefore, the average field at 
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Too = 

(R)=^T(r,-r,).(m,), (10) 

j 

(T is the dipolar tensor) is expected to approach zero. 

One needs to mention tliat the orientational structure 
of surface water around spherical cavities and non-polar 
solutes results in a positive value of the cavity electro- 
static potential] ^^'^^ The potential is apparently rather 
uniform resulting in (R) ~ in both the present simula- 
tions with the off-center dipole, as well as in our previous 
study of the same Kihara solute with the dipole placed 
at the solute's center 

In contrast to the first-order orientational order pa- 
rameter p\ the second-order orientational parameter 

pi = {2N'r' J2 [3(m,-f,)^-l] (11) 

j=l,Ni 

is of the order — —0.19, indicating that about 40% of 
first-shell waters have their dipoles directed perpendic- 
ular to the surface normalJ^"— This parameter decays 
only slightly upon the insertion of the solute dipole (Fig. 
II. 



0.04 
0.036 




RDIO 
.. .RD9 


0.032 






2 4 6 8 10 




1^' 


' ° , ■■■■■■■■ 
■■ 





FIG. 4. Second [pi, Eq. ^] ( main panel) and first [pi, 
Eq. ([9])] (inset) orientational order parameters for RD9 (r^ — 
9 A, filled squares) and RDIO (r<i = 10 A, open squares) 
configurations. 

Taken together, the results for p{ 2 suggest that there 
is little asymmetry in the distribution of the interfacial 
dipoles. The condition (i?) ~ in Eq. (fTO|) then implies 
symmetric J-{0, R), which can be achieved with a reduced 
number of model parameters: R2 = — i?i and Fi = F2 ~ 
Fq. The free energies of the three states then become 

Fg{R) =RV{2k,) 

F^{R)={R-R^f/{2K^) + Fo (12) 
F2{R) ={R + Ri)y{2Ki) + Fo 

Four model parameters Kg, ki, Ri, and Fq were fitted 
to (t|.(too) from MD data. They are listed in Table |T] 
and are used to produce plots of i?o(TOo), cr^(TOo), and 
J-(mQ, R) presented in the Results section below. 



TABLE I. List of model parmeters produced by fitting the 
three-state model to aj^{mo) from MD simulations. = 
{koi)~^^^, a = g, 1 is the characteristic length of penetration 
of water's surface fluctuations into the solute [Eq. ((6|]. 

Configuration rg, A pFp t% k l^i^, k Ri, V / k 
RD9 9 7l 3l) O O 

RDIO 10 7.2 1.5 1.2 1.0 

III. RESULTS 

All the simulation results that we discuss below have 
been produced for two distances between the point dipole 
in the solute and its center (Fig.H]), = 9 A and Vd ~ 10 
A. We will reference them throughout below as states 
RD9 and RDIO, respectively. 

There is a fairly broad intermediate region between 
two linear scaling regimes of the reaction field i?o with 
the solute dipole toq. This intermediate, approximately 
quadratic, scaling is a signature of transitions between 
the low and high free energy states. This is clearly seen 
in the progression of free energies T(rno,R) with chang- 
ing Too shown in Fig. [S] for the RDIO state. The solid 
lines in the figure are obtained directly from distribu- 
tions of R from MD simulations, while the dashed lines 
are the fits to the three-state model outlined above. The 
overall picture of the alteration of the free energy surface 
with the external dipole is remarkably consistent with 
the sketch shown in Fig. [T)d pointing to a discontinuous, 
first-order character of the cross-over between the ground 
and excited surface states. 

As is expected from the general phenomenology of dis- 
continuous transitions , ^1^"^'^^ the fluctuations of the or- 
der parameter are non-Gaussian close to the transition 
point, reflecting transitions of the system between the 
two alternative states. Such enhanced fluctuations are 
typically reflected by spikes of variances, such as the heat 
capacityi^^ Our observations are no different and the vari- 
ance cr^ = (((5i?)^) shows a peak at the transition value 
of the solute dipole moment (Fig. [6)d) . The results are 
plotted against a dimensionless dipole magnitude that re- 
flects the ratio of the characteristic held produced by the 
solute to the one existing within the solvent 

Too = — , (13) 

rris \rdsj 

where rUs and as are the solvent (water) dipole moment 
and diameter, respectively. In addition, rds = — rns 
(Fig. [2|) is the distance between the point dipole and the 
solute- water interface. Similar ideas are used to intro- 
duce the reduced instantaneous held 

R* - Rrl/m, (14) 

the average, Rq, and variance, {{SR*)'^), of which are 
shown in Fig. [S] 

The spike in the fleld variance is accompanied 
by a strong increase in the relaxation time of the 



4 




FIG. 5. Free energies T{mo,R) [Eq. for the RDIO state 
(rd = 10 A) at mo values indicated in the plots. The reac- 
tion field in the abscissa is in V/A; the dashed lines refer to 
the fits to the three-state model [Eqs. ([12}]. The model 
parameters used in the plot are listed in Table |T] 



normalized time self-correlation function, S(t) = 
{5R{t)5R{0))/{{SRf). The relaxation time shown in Fig. 
[6}: was calculated as Toff = dtS{t). The appearance of 
the peak in the relaxation time is reminiscent of critical 
slowing down well established for phase transitionsr^i^ 
but how far that language can be carried over to our 
situation is not entirely clear. 

The preceding discussion clearly indicates that fluc- 
tuations of water's electric field are expected to be non- 
Gaussian in the range of solute dipoles allowing switching 
of surface waters between the ground and excited states 
in the spirit of Figs. [1] and El We indeed observe a spike 
in the nonlinearity parameter xg [Eq. in the tran- 
sition region (Fig. [7^). The extent and the breadth of 
this non-Gaussian behavior are additionally illustrated 
in Fig. [7b where the direct calculation of the stiffness pa- 
rameter K from the reaction field, k = Ro/niQ, (circles) 
is compared to the expectation from the linear response 
according to Eq. k — /3(t^ (squares). Since the linear 
response function is independent of the dipole moment, 
the horizontal dashed line in Fig.[7|D indicates the value to 
which both calculations would converge if linear response 
held. In order to provide an additional reference point to 
our results, the horizontal dash-dotted line shows k for 
a point dipole placed inside a cavity in dielectric contin- 
uum. The reaction field in this approximation is derived 
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FIG. 6. The average, Rq, of the reduced field R* as defined 
by Eq. (|14|) (a) and its variance (b) vs the reduced solute 
dipole defined by Eq. (|13|) . The filled squares and solid lines 
in (a) and (b) refer, respectively, to the MD results and their 
fit to the three-state model for the RD9 configuration; the 
open squares and the dashed lines carry the same information 
for the RDIO configuration. Panel (c) shows the relaxation 
time Tcff of the water field self-correlation function (see text) 
obtained from MD simulations for RD9 (solid line) and RDIO 
(dashed line) configurations. 



in the SM2i 

Ro = 



Too e — 1 
^2e + 1 



1 



(15) 



where z = rd/ros- The continuum cavity radius in this 
calculation was chosen at ro^ = 11.5 A, slightly below 
the position of the first peak of the solute-solvent pair 
distribution function at ~ 12 A(Fig. [5]). 

The results presented here pose a significant question 
of the nature of the excited state of the surface water. A 
related question is what is the spatial extent of the ex- 
cited water domain. The results for the alteration of the 
orientational order parameters with changing toq ( Fig. 
HI support the notion of a relatively localized perturba- 
tion, not spreading into the entire hydration layer. 

In order to address the question of the water struc- 
ture in a domain nearest to the solute dipole we have 
separated a solid angle of 27.4° drawn with the dipole 
direction as the axis and cutting a portion of the surface 
circumference equal to 2a s (see SM^i for a drawing) . The 
pair solute- water distribution function was calculated for 
waters residing within this solid angle only and is shown 
in Fig. [51 A new peak emerges from the nearest solute- 
water distance with increasing toq. The area under the 
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FIG. 7. Nonlinearity parameter XG [Eq. (HJ] (a) obtained 
from MD simulations (points) and from fitting the three-state 
model (lines). Panel (b) shows the reduced stiffness parameter 
i^'f'os = {ros/lnf' calculated in different approximations for 
RD9. The dash-dotted line, k = Ro/mo, is from the dielectric 
continuum calculation for the reaction field at a dipole inside 
a dielectric cavity [Eq. (|15p ]. The dashed line shows k — I3a\ 
obtained from the variance a\ at mo = 0. The open circles 
show « = Ro/mo and open squares refer to k = /3cr|j, both 
from MD simulations. In linear response, circles and squares 
are expected to collapse on the dashed line, Eq. @. 
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FIG. 8. Solute-oxygen pair correlation function for RD9, in- 
side the solid angle defined in the text (see also the SM^^-). 
The numbers in the plot indicate values of the dipole moment 
mo. The curves are shifted vertically by 1.0 for a better view. 



peak shows that exactly one water molecule breaks from 
the surface hydrogen-bond network, shifting toward the 
solute. This interpretation is consistent with the exci- 
tation free energy of Fq ~ 7 k^T (Table close to the 
enthalpy required to break two hydrogen bonds in SPC/E 
water 

The creation of the excitation defect in the hydration 
shell results in the frustration of the nearest bonds net- 
work seen in the alteration of the distribution of 0-H-O 
angles of first-shell waters within the same solid angle as 
used in producing Fig. [S] These distributions at chang- 
ing mo are shown in Fig. [9] Similar to the previous ob- 
servations by Sharp and co-workers,— the population of 
the buckled, large-angle state grows with increasing mp. 
This is a behavior generic for polar (ionic or dipolar) so- 
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FIG. 9. Distributions of O-H-O angles 6 for the RDIO con- 
figuration at mo equal to 0.01, 3.5, and 10 D. The distribu- 
tions were calculated according to the algorithm suggested by 
Sharp et ai^ in the region of the first hydration layer within 
the solid angle cutting the circumference length of 2(Ts at the 
surface (see also the SM^^). The inset shows the definition 
of the bond angle and the distance cutoff within which the 
angles were sampled. 



lutes frustrating water's hydrogen-bond network by their 
electric fields. 



IV. DISCUSSION 

We have studied the statistics of electric field fluctua- 
tions produced by surface water inside the solute. The 
main finding of this study is the realization that the 
statistics can be profoundly altered by placing a dipole 
next to the interface. With increasing dipole magnitude, 
waters closest to the dipole alter their hydrogen-bond 
network in a way that creates a domain of surface waters 
in an excited state, with the free energy lifted by approx- 
imately the energy of breaking two hydrogen bonds (~ 7 
/cbT, Table m. A crude estimate of what is needed to 
shift the system into this excited state in provided by 
Fig. [61d: a dipole twice as large as water's dipole placed 
one water diameter from the interface puts the system 
at the cross-over point at the peak of non-Gaussian fiuc- 
tuations. This outcome seems to be quite generic: once 
surface particles have close probabilities to occupy alter- 
native states, the variance of the fiuctuations is in excess 
to the expectation from the linear responsci^ 

The local nature of the excitation, most likely limited 
to one water breaking from the network and creating a 
corresponding elastic deformation around it, contributes 
to a significant range of dipole moment magnitudes at 
which the system finds itself distributed between two 
local minima, each characterized by a linear response. 
This property of the local surface transition, clearly dis- 
tinct from much more localized, in the parameters space, 
phase transitions^ will allow one to observe the effects 
discussed here for a substantial range of solute perturba- 
tions. 

The surface defect and the corresponding frustration 
of the network characterizing the excitation do not nec- 
essarily require a surface dipole. They might be achieved 
by other causes such as ionic or hydrogen-bonding sur- 
face sites. An optical probe placed near such a domain 
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of waters in the excited state will record an optical re- 
sponse distinct in both the spectral shift and the width 
from probes next to the ground-state interface. In this 
regard, the inhomogeneous spectral broadening, which is 
admittedly hard to extract experimentally, is neverthe- 
less a more sensitive probe of the local interfacial struc- 
ture than the spectral shift. The Stokes shift dynamics 
can also slow down in the transition region, as is seen in 
Fig. A mosaic of ground- and excited-states domains 
at an extended interface of a nanometer-scale solute will 
then contribute to energetic and dynamic heterogeneity 
of the electric field fluctuations However, the dynam- 
ics of the electric field observed here (ps) are still sig- 
nificantly faster than that found for hydrated proteins 
(ns) .— The three orders of magnitude difference sug- 
gests a possibility of cooperativity between many surface 
excitations at proteins' polar/ionic surface sites. Alterna- 
tively, the excited surface state of a protein might still be 
localized, and the slow dynamics come as a result of pro- 
tein motions modulating the interfacial polarization)^ 
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